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SUMMARY 

Rotor-support systems interaction with parametric excitation is considered 
for both unequal principal shaft stiffnesses (generators) and off-set disc 
rotors (ventilators). Instability regions and types of instability are com- 
puted in the first case, and parametric resonances in the second case. Com- 
puted and experimental results are compared for laboratory machine models. 

A field case study of parametric vibrations in industrial ventilators is re- 
ported. Computed parametric resonances are confirmed in field measurements 
and some industrial failures are explained. Also the dynamic influence and 
gyroscopic effect of supporting structures are shown and computed. 

INTRODUCTION 

With growing rotor speeds or machine capacities the influence of support- 
ing structures on stability and critical speeds is of increasing importance, 
specially when some rotor parameters are variable. This is often the case in 
machines used in mining and power generating industries, for example, large 
generators, industrial ventilators, and turbopumps. The unequal rotor princi- 
pal stiffnesses in the first instance and unequal principal moments of inertia 
of the ventilator disc in the second instance, with the flexible anisotropic 
support interaction, cause the parametric effect* 

^ 120-MW generator of Do lmel- Wroclaw production is shown in figure 1. 

Its rotor is of 9.157 m length, with a span of 6.858 m between bearings. The 
wiring slots distribution gives rise to small unequality of the rotor principal 
stiffnesses k x ,ky. Two types of mine ventilators are shown in figure 2: 
large ventilator No. I of the type SLM (fig. 2(a)) of nominal output of 
320 m/min, and universal ventilator No. II of type FKD-30 (fig. 2(b)). The 
first has an overhang shaft of full length 0.908 m. The disc is between bear- 
ings in the second ventilator with a shaft of 0.773 m length. Production errors 
cause the unequality of principal moments of inertia of the disc. 

GENERAL THEORY 

The rotor-support system can be considered as two connected and dynamically 
interacting subsystems. One (rotor) is a discrete parametric system (fig. 3, 
bold lines). The second is nonparametric, contains the whole supporting struc- 
ture, and is usually very complicated. A receptance (impedance) matrix is hence 
used for description of its dynamic behavior. Its elements can then be mea- 
sured if not calculated. All (usually infinity) but the connecting coordinates 
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of the supporting structure can be eliminated (ref. 6) from the equations of 
motion of the systems, giving their general matrix form 



1 

P 



X= G + Psin ot ♦ Q cos <yt 

H K H K 


( 1 ) 


Here all the matrices are of the order equal to the sum of the number of de- 
grees of freedom of the parametric subsystem (rotor) and the number of connect- 
ing coordinates. In the cases considered, for instance, the last number is 
four (fig. 2); y r and y u are, respectively, real and imaginary parts of the 
dynamic impedances of the supporting structure along the connecting coordinates; 

p j and Q are force vectors. For an offset disc mounted flexibly (fig. A) 
or for a rotor of different principal stiffnesses (the two cases are here con- 
nected only formally) , the equations are 

For a rotor of different principal stiffnesses 

M(t)=M; C(t)*C; Klt)*K>£(Ecos‘frt*Fsin #t) (2a) 


For an offset disc and circular shaft 
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where S^, — are the first and second moments of the rotor mass with re 
spect to the coordinate planes connected with the rotor; D r , D n , D r are 
corresponding deviation moments; and <= n t 


b=!(b„*-b ?5 ) 


Finally 
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where p and 0 are the arguments of the solutions 


x = e iJUt 2(a, sin -jt + b, cos ^ t) 


( 10 ) 


of the homogenous part of equation (1) at stability limits. 


Instability 

The conditions of existence of steady solutions (10) , giving the boundary 
equations in the stability region, in the equivalent forms for odd or even 
numerals i, are 



= o 


do 


Consecutive approximations are shown with broken lines in equation (11). The 
first approximation is 
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= O ( 12 ) 



The instabilities occur in the vicinity of circular frequencies 

X . + X. 

e., — * 

jk n 


n - 1 , 2 , 3 , . . . 


(13) 


The width of the instability region depends on e and damping C. Here Aj 
are the natural frequencies of the system, being the roots of equation (12) at 
zero parametric excitation (e = 0) and no damping (C — 0, Y u = 0) . 

Resonance Speeds 

Forced vibrations analysis gives for the resonance speeds (refs. 3 

and 6) 


and 


(2n - l)« cr = 


2n ii) = X (oi ) 
cr z cr 


~\ 


(14) 


Here X z is the axial natural circular frequency of the system. 

LABORATORY MODEL TESTS AND COMPUTATIONS 


Computations were done for a laboratory model (fig. 5) and compared with 
the tests results. 


Generator Model - Instabilities 

The instability regions were computed for the model with the shaft of dif- 
ferent principal stiffnesses k^ky and rotor block between bearing (no disc 5 
in fig. 5). The determinants (12) give equations of instability region bound- 
aries in the form |w(y,9)| = 0. Their roots Ui,6i define the stability 
limits. An example of computed results is given in figure 6. This diagram 
corresponds to the second-order combined instability region 0.5(9 X ^ + 9y]_) 
around point N (fig. 6). The solutions of the equation 

iw O,o)l = o 
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for zero parametric excitation (e = 0) and no damping (C = 0 y =0) are 
straight lines ' u ' 


p (± L if -Lis = x 

r 2 j 


intersecting at point N (fig. 6). For finite parametric excitation the lines 
curve and do not intersect, separating in the direction 0 for frequencies 
sums l/n(Aj + A k ) or in the direction y for their differences l/n(Ai - Av). 
Hence there are no instabilities (also with no damping) corresponding to the 
frequencjr differences. The curves in figure 6 correspond to constant values of 
the determinant, and its zero values (denoted with the crosses) give the in- 
S ^ a region. They approach each other with increasing damping and, over 

1 K S r ll '? ltln ? V f lue ’ instability disappears. This was traced with changing 
s ait length AL, and computed and experimental results are compared in fig- 
ure 7. Corresponding oscillograms for the x and y directions are given in 
igure 8. The instability region discussed here is denoted by 1 in the insta- 
liity regions diagram (fig. 9). There are computed instability regions (with 
e type of instability) shown in figure 9 and against them the observed experi- 
mental instabilities shown with the crosses. 


Ventilator Model - Resonance Speeds 

Analysis of forced vibrations was done for the model with circular shaft 
and the disc 5 (fig. 5) slightly offset. Hence B n? and are not equal 

ana D ? as well as disbalance forces P and Q are not zero, but 

e 0. Graphical solution of equation (1) for the laboratory model is given in 
figure 10. The continuous curved lines are the computed natural frequency 
curves. Abscissas of their points of intersection with the straight lines 

w (principal), A = 2 oj, A = 3to . . . give the parametric resonance speeds, with 
one, two, three, ... vibration cycles per one revolution. Computed values of 
10 resonance speeds (taken from diagram, fig. 10) in s~l are 


“crl 20 - “d2 " 36 ' "cr3 ‘ 46 ' “ c r4 ' 58 , <0^ . 82, 


“cr6 * ”• “cr7 “ 108 > “ cr8 ’ 112 > “ cr9 ' 140. » crl0 > 175 


Experimental resonance diagrams measured for the right support (fig. 5) are 
given in figure 11 as continuous broken and dotted lines for the directions 
y, x, z, respectively. They compare well with the computed values from fig- 
ure 10. The vibrations recorded for speeds close to resonance, numbered in the 
diagrams (figs. 10 and 11) are shown in figure 12. It can be seen that close 
to resonance speed of the type no) cr = A i (m cr ) vibrations of frequency nw domi- 
nate. Comparison of figures 10, 11, and 12 shows good consistency of theoreti- 
cal and experimental results. 
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FIELD CASE STUDY 


The above results were applied in field studies of vibrations of industrial 
ventilators that developed repeated failures (shaft breakages). These were 
ventilators No. II of the type FKD-30 (fig. 2(b)) mentioned in the INTRODUCTION. 


Analysis 

The model of the ventilator consists of two subsystems: the discretized 

rotor (fig. 13) described in table 1 and the supporting structure containing 
the ventilator body with bearings and foundation (fig. 14) . The supporting 
structure, for its dynamic interaction with the rotor, was dynamically defined 
with its receptance matrix along the connecting coordinates x^, yp, x£, y2 


xlxl 

C x1x2 

C x1y1 

c x1y2 

x2x1 

C x2x2 

c x2y1 

c x2y2 

ylxl 

C y1x2 

c yiyi 

c yiy2 

'y2x1 

C y2x2 

c y2y1 

c y2y2 


(15) 


The receptances c ik were measured on site. The measuring scheme is given in 
figure 15. The notations in figure 15 are 1 - ventilator type II, 2 - electro- 
dynamic exciter, 3 - displacement pickup, 4 - generator, 5 - vibration measuring 
apparatus, 6 - two-ray cathode oscilloscope, 7 - frequency meter, 8 - phase 
meter, 9 - ammeter. 


The exciting force was measured indirectly (by current measurement). Re- 
ceptance diagrams are given in figure 16. Because of symmetry all the cross 
receptances for the perpendicular directions are zero. The other equal pairs 

„„„ _ „ r . . c .. where i = 1, 2. On the basis of the re- 

are denoted c iiy , c ij x » c ijy’ wueJ - e ’ 

ceptance matrix (15) the impedance (dynamic stiffness) matrix 


[Y ik ] = [G ik ] 


-1 


was calculated and introducted in the computing program. Forced vibrations 
were found as partial solutions of equation (1) for the case of equal rotor 
stiffnesses by application of the small-parameter method. Their analysis gives 
the formulas (14) for the resonance speeds. Graphical solution is given in 
figure 17. The frequency lines are the solutions of the homogenous part of 
equation (1). With the approximation used, there are computed two principal 
resonance speeds (in s"l) , 344 and 280; two subharmonic ones of second order 
(of two cycles of vibration per revolution), 172 and 142; and two of third 
order, 114 and 94. 
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Experimental Measurement on Site 


^ perimental resonance test was made on site (in the mine) in the ven- 

ure lS 7 * ? COndltion * . The measuring arrangement is shown in fig- 

1 denotes the ventilator tested; 2 denotes the dc motor; 3 de- 

inrapiaLLs nS 5°^e r tr ^ sduce ^; 4 denotes the vibration measur- 

ing apparatus, 5 the cathode ray oscilloscope; and 6 the frequency meter. 

The resonance diagram for the supports for the practically allowed speed 
range from approximately 1200 to 2100 rpm is given in figure 19. The two ob- 
served subharmonic resonances, 138 and 174 8 -l, are in this speed rangelboth 

172 a"?? ° are Cl ° Se t0 the corres P° ndin 8 computed values, 142 and 


CONCLUSIONS 

n , Tha second -order subharmonic of frequency 142 s"l coincided with the nomi- 
nai working speed of ventilator II, which was below 1500 rpm. Hence it was the 
cond-order parametric resonance that was responsible for the shaft breakages, 
upporting structure appreciably lowers critical speeds. The gyroscopic effect 

"th a Le 0 COnS L ered and COmPUted ‘ Tt 15 SpeCi3lly Visible in ^ ventilators 
d ha ° varha ng rotor. The frequency diagram for that type of ventilator, 
gnated No. I in the INTRODUCTION, is given in figure 20. Important 

! ffGCt lnf luence is visible. It was the principal resonance at 

that w ^ s res P° nsible for the failures in those ventilators with nominal 
working speed below 1500 rpm. 
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Pig. 5 Scheme of model stand 
1-DC motor, 2-elastic coupling, 

3 -shaft ,4-rotor mass, 5 -disc, 6 -hall 
hearing, 7-anisotropic he coring, 
8-flat springs, 9-inductive sensor, 
10-piezo-electric sensors, 1 1 -re- 
corder, 12-recorder tape,13-ele- 
ctrodynamic exciter, 14-generator 
of harmonic vibration, 15-cathode 
oscilloscope # 



Pig. 6 
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Fig. 18 
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